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Army. 


I v*T Where copyrighted material is quoted, permission has been 
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distribution is quoted, permission has been obtained to use the 
material. 


j (7( s Citations of commercial organizations and trade names in 
this report do not constitute an official Department of Army 
endorsement or approval of the products or services of these 
organizations. 

A iL In conducting research using animals, the investigator(s) 
adhered to the "Guide for the Care and Use of Laboratory 
Animals," prepared by the Committee on Care and use of Laboratory 
Animals of the Institute of Laboratory Resources, national 
Research Council (NIH Publication No. 86-23, Revised 1985). 


f u(_^ For the protection of human subjects, the investigator (s) 
adhered to policies of applicable Federal Law 45 CFR 46. 

In conducting research utilizing recombinant DNA technology, 
the investigator(s) adhered to current guidelines promulgated by 
the National Institutes of Health. 

aM- In the conduct of research utilizing recombinant DNA, the 
investigator(s) adhered to the NIH Guidelines for Research 
Involving Recombinant DNA Molecules. 

In the conduct of research involving hazardous organisms, 
the investigator(s) adhered to the CDC-NIH Guide for Biosafety in 
Microbiological and Biomedical Laboratories. 
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i INTRODUCTION 

Gadolinium-based magnetic resonance (MR) contrast agents offers exceptionally high 
sensitivity for the detection of breast cancers [1-6]. The specificity of gadolinium- 
enhanced MR mammography to distinguish benign from malignant tissues, however, has 
been limited by the fact that some normal and benign masses also exhibit substantial 
contrast enhancement [3,5,6]. Alternatively, properties related to the "temporal rate" of 
contrast uptake has being studied as a potential complimentary discriminator of benign 
from malignant breast neoplasm [7-11]. Investigators have noted, malignancies tend to 
exhibit a rapid contrast enhancement profile that is presumably due to an anomalous 
vascular supply and highly permeable capillaries. Non-invasive methods to detect and/or 
discriminate between normal, benign, pre-malignant and malignant tissues based on 
vascular alterations would have significant diagnostic value, particularly in cases that are 
a challenge for conventional mammography such as in dense breasts. 

Acquisition and processing methods to quantify MRI contrast agent enhancement kinetics 
have been widely variable across studies [7-11]. For accurate characterization of contrast 
uptake features of breast lesions in a clinical setting, critical parameters include: temporal 
resolution, spatial resolution, volume of tissue coverage, and the method to quantify 
signal changes. Information regarding “rate of enhancement” is obscured if the temporal 
sampling is inadequate to record details of the lesion enhancement curve such as the point 
contrast onset and maximal slope. More serious errors occur if the imaged volume only 
partially includes an enhancing lesion or does not include all lesions in the breast. These 
issues arise, in part, due to the well known trade off between volume coverage, spatial 
resolution, and temporal resolution. 

The central objectives of this project are technical refinement and clinical 
implementation of new dynamic contrast-enhanced MRI methods designed to provide 
full volume coverage of both breasts at good spatial and temporal resolution for detailed 
quantification of dynamic contrast enhancement properties [10]. Key elements of this 
approach involved specialized software for data acquisition (i.e. MRI pulse sequence), 
motion correction and image reconstruction, and kinetic analysis of contrast-enhanced 
tissues. The methodology developed for this project allowed dynamic scans through 32 
sections of both breasts acquired at 10-13 second temporal resolution - this was a 4 fold 
higher temporal sampling than the commercially available sequence. Additional post¬ 
processing, referred to as viewsharing, further reduced temporal sampling to 5-6seconds. 
This technique was applied to study 102 women in whom a breast abnormality was 
detected by mammography, physical exam, or ultrasound. 
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METHODS 

MRL Data Acquisition 

Breast MR exams were performed on 1.5Tesla General Electric Signa systems using a 
commercial (GE) bilateral phased-array breast coil with insert pads for mild 
anterio/posterior (A/P) compression. The MRI protocol for all patients was as follows: 
Series 1: Axial Locator 

Axial sections 10mm thick; 2D Spin-Echo (SE); Echo-Time (TE)=12ms; Rep-Time 
(TR)=300ms; Matrix = 256x128; Field-of-View (FOV)=35cm; Excitations 
(NEX)=1; ScanTime = lmin 29sec. 

Series 2: Axial T2-weighted 

Axial sections 6mm thick; 2D Fast-Spin-Echo (FSE, 8echo-train); TE=102ms; 
TR=3000ms; Matrix = 256x192; FOV=28-32cm; NEX=2; Fat signal suppressed 
(FatSat); Freq axis=Ant/Post; ScanTime = 5min 23sec. 

Series 3: Axial T1-weighted Pre-Contrast 

Axial sections 6mm thick; 2D; Phase-Offset-Multi-Planar; TE=14ms; TR=400ms; 
Matrix = 256x192; FOV=28-32cm; NEX=2; Fat signal suppressed (FatSat); Freq 
axis=Ant/Post; ScanTime = 5min 20sec. 

Series 4A: Coronal 3D T1-weighted Pre-Contrast 

Coronal sections 2-5mm thick; 32slices; 3D Spoiled Gradient-Recalled Echo 
(SPGR); TE=3.3ms; TR=10ms; Flip=40°; Matrix = 256x128; FOV=28-32cm Rt/Lt, 
FOV=14-16cm Sup/Inf; NEX=4; ;Freq axis=Rt/Lt; ScanTime = 2min 38sec. 

Series 4B: Dynamic Coronal 3D Tl-weighted 

Repeat 4A with Matrix = 256x32 and NEX=1 in each pass; 20 passes, then one pass 
of Matrix = 256x128; Scan Time = 4min 6sec. 

Series 4C: Coronal 3D Tl-weighted Post-Contrast 
Repeat 4A Scan Time = 2min 38sec. 

Series 5: Axial Tl-weighted Post-Contrast 
Repeat Series 3; ScanTime = 5min 20sec. 

The dynamic series was built upon a 3D-volumetric rf-spoiled gradient recalled echo (3D 
SPGR) sequence since it (a) offers excellent sensitivity to gadolinium-induced changes, 
(b) produces high signal-to-noise images, and (c) provides full volume coverage of both 
breasts with contiguous image sections. Series 4 segments A, B, and C were 
concatenated for speed and to eliminate need for machine adjustment; this helped 
minimized patient movement. A standard dose of gadolinium-contrast material 
(0. ImMol/kg; Magnevist, Berlex Labs, Wayne, NJ.) was manually injected as a bolus at a 
specified time during segment B. 

MRI Data Processing 

Scan time and spatial resolution are proportional to acquisition matrix size. Dynamic 
changes (due to gadolinium), however, are most heavily concentrated near the center of 
the 3D Fourier-encoded MRI dataset referred to as “k-space”. Therefore, to successfully 
record most dynamic contrast changes, while greatly reducing scan times for high 
temporal resolution (i.e. a 3D image set every 10 seconds), only the central k-space data 
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were acquired during the gadolinium injection (i.e. Series 4 B). The lost spatial 
resolution was partially restored by reuse of previously acquired peripheral k-space data 
(i.e. from Series 4 A). This process of splicing high temporal resolution central k-space 
data blocks with high spatial resolution peripheral k-space data is referred to as the 
“keyhole” method, and was extended to 3D specifically for this study [10]. 

An additional a post-processing technique referred to as “viewsharing”, was applied to 
effectively double the temporal sampling rate to a 3D image set every 5 seconds. 
Viewsharing is the synthesis of data blocks at each intermediate timepoint by 
combination of the latter half of a given data keyhole block with the initial half of the 
subsequent data block. As with each original keyhole data block, synthetic viewshared 
blocks were submitted for keyhole reconstruction. In this study, a series of 21 keyhole 
blocks were acquired and interleaved with 21 viewshared blocks for a total 42-point 
timeseries spanning 4-5minutes. 

Contrast enhancement analysis requires that the patient remain relatively still during all 
of Series 4; often this was not the case. A simple linear motion correction algorithm was 
incorporated within keyhole/viewshared image reconstruction to reduce motion artifacts. 
The algorithm is described in detail elsewhere (Medical Physics, in press; see appendix), 
but briefly it involves matching the phase of each k-space block to the pre-contrast k- 
space data. In effect, this phase adjustment is equivalent to a 3-dimensional translation of 
each imageset toward spatial coregistraion with pre-contrast images. Since the pre¬ 
contrast data also serves as a subtraction mask, this correction reduced motion artifacts 
that are most apparent on subtracted images. Motion-corrected, viewshared, keyhole 
anatomic and subtracted images were viewed cinegraphically on the SUN SPARC 10 
workstation where all reconstruction/analysis was performed. Rapidly enhancing lesions 
are most conspicuous on the cinegraphic subtracted images. This CINE format aided the 
user to electronically define a region of interest (ROI). 

Lesion ROI enhancement vs time curves were automatically fit to exponential saturation, 
and exponential saturation-decay models. These mathematical functions are derived from 
compartmental models that describe the temporal distribution of contrast media between 
tissue vascular and extracellular spaces. These models, albeit simplistic, suggest that 
capillary permeability and surface area (P»S product) are most relevant to early 
enhancement rate. Enhancement rate is defined as [fitted enhancement time, tau]’ 1 . Also, 
overall enhancement amplitude (i.e. degree of signal change) reflects size of extracellular 
space and amount of contrast leakage. Fitted enhancement amplitude and temporal 
properties were recorded on disk for each ROI. Image excerpts, contrast enhancement 
curve, and fitted results for an invasive carcinoma are illustrated in Figure 1. 
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Invasive Ductal Ca 



Figure 1. Excerpts from an enhancement time series of a invasive ductal carcinoma in the 
upper-outer quadrant of the right breast. This lesion exhibits a high enhancement rate 
(P»S ~ l/18sec = 57xl0' 3 sec') and washout; both signs are strongly suggestive of cancer. 

Patient Recruitment - Image Review 

Recruited patients were referred from the Breast Care Center of the University of 
Michigan Comprehensive Cancer Center. Entrance criteria includes any patient with a 
breast abnormality detected by mammography, physical examination, or ultrasound, and 
which will be biopsied. Prior to MRI study, all patients had a routine clinical evaluation 
which included mammography/ultrasound and breast physical examination. 

Images/films/reports were reviewed by radiologists (project co-investigators) in the order 
of: (a) conventional mammography images (Xray, ultrasound) and previous reports; (b) 
conventional MRI including pre- and post-gadolinium; and (c) subtracted MR images 
(pre-gad subtracted from post-gad). The radiologist defined which slice/region best 
corresponded to the lesion/tissue sampled by surgical/core biopsy. Occasionally, 
additional focal tissues were noted on MRI that were correlated with mammographically- 
stable (>lyear) lesions. These were also included in the analysis as benign lesions. 
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Statement of Work Tasks and Status 

SOW Task 1: Patient recruitment - uniform over the four years of the study 
One-hundred breast patients will be studied during this four-year prospective 
imaging trial. Volunteer patients will be recruited from the Breast Imaging 
Section of the Breast Care Center at the University of Michigan Hospital. 
Entrance criteria will include patients with a breast abnormality detected by 
mammography, physical exam, or ultrasound which will be biopsied. 
Malignancy, grade, type and lymph node involvement will be determined 
for each patient by pathologic reading. Approximately 30 malignant and 70 
benign lesions are expected from this population. 

Status Task 1: Originally, one-hundred women were to be recruited for this study. 
Ninety-four women were enrolled and scanned during project period; two failed to 
complete the MRI exam due to claustrophobia. An additional ten women were included 
from the pilot study completed prior to the project period. The same scan/processing 
techniques applied and these women had biopsies. This provided data from a total of 102 
subjects. Per the stated selection criteria, any patient with an abnormality sufficient to 
warrant biopsy was considered for this study. This was done to not limit selection to only 
those with highly suspicious discrete masses. As such, a relatively large number of 
patients were included in the study by being referred for architectural distortion and/or 
calcification clusters, but without having a discrete mass. In some of these cases focal 
enhancement in the region of mammographic concern made correlation of biopsy 
sampling and MRI ROI definition straightforward. In cases lacking focal enhancement in 
the region of mammographic abnormality the correlation of the biopsy/surgical site on 
the MRI was more difficult. In 27 cases, the radiologist used the mammogram and MR 
images to provided a "best estimate" of the biopsy area on the MRI even though the MRI 
exhibited non-focal or moderate parenchymal enhancement. All of these 27 cases would 
have been read as "negative" by MRI due to lack of focal enhancement. The pathology 
results on these 27 cases were all non-malignant - most common path reading was 
fibrocystic change, with or without calcifications. In the interest of testing the value of 
the quantitative MRI parameters, data from these ROI's were maintained in the statistical 
analysis described in the Results. At the project end, quantitative enhancement data were 
collected from total of 39 cancer and 54 non-cancer tissues. 


SOW Task 2: Enhancement rate-constant determination - each patient 
Each patient will be studied by a new dynamic MRI sequence during bolus 
administration of gadolinium-DTPA contrast. The data acquisition and 
processing sequence will generate a temporal series of 32 data sets, where 
each set represents 32 tomographic slices through both breasts, acquired at a 
rate of 1 set per 10-13 seconds. Existing specialized interactive software 
will be employed to retrospectively view each tomographic slice 
cinegraphically to identify gadolinium enhanced tissues. Enhancement rate- 
constants will be quantified by numerical fit to signal change versus time 
using a exponential saturation model. 
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Status Task 2: Complete. 

SOW Task 3: Software development: View Sharing & Automated Fitting - 1st year 
View sharing is a post-processing concept designed to effectively double the 
dynamic sampling rate using existing data. Software to execute and test this 
concept will be written for simulated and actual MRI data. Pending 
satisfactory performance of view sharing, archived data sets will be 
reprocessed, and subsequent data sets will be processed with view sharing 
thus yielding 64 data sets at 6sec dynamic temporal resolution. Addition 
development work will be directed toward export of datasets from 
workstation back into MRI system database for more portable viewing. 

Status Task 3: Complete. Keyhole reconstruction, viewsharing, and motion-correction 
algorithms were applied to each case. Total number of viewshared datasets, however, 
was limited by computer memory to 42 timepoints instead of 64. We do not believe this 
change affected results since the critical "rapid-enhancement" segment of the curve 
happens early in the timeseries. We elected to retain dynamic breast image datasets on 
the workstation since transferal and storage of large image sets (>2000 images per case) 
is cumbersome. Instead, we added a high-quality hardcopy unit and a 2.6Gbyte magneto¬ 
optical drive to the SPARC 10 workstation, thus harcopy/analysis on the MRI system was 
not essential. 


SOW Task 4: Statistical summary - final two months of study 

Quality of curve fitting and appropriateness of selected model will be statistically 
assessed on individual patient data. This will be performed throughout the study, as 
will assembly of other diagnoses and pathological outcomes into a database. Final 
inter-patient statistical summaries will occur in the last few months of the study. 
Quantitative enhancement rate constants, conventional pre-biopsy diagnosis, and 
pathologic characterization will be submitted for statistical testing of stated 
hypotheses regarding diagnostic sensitivity, specificity, positive and negative 
predictive values of enhancement rate constants relative to conventional pre-biopsy 
diagnosis. Classification power for tissue group type, and correlation with tumor 
grade, size, invasiveness, and lymph node involvement will be determined. 

Status Task 4: Near Complete. Statistical analyses are complete for MRI parameter fits 
to: (a) enhancement time, t, (time-constant of the exponential enhancement curve); (b) 
enhancement rate (1/t) related to permeability surface-area product; (c) enhancement 
amplitude; and (d) product of enhancement rate and amplitude for available ROIs. Since 
these are objective numerical classifications of cancer vs non-cancer, a Receiver 
Operating Characteristic curve was developed for each MRI parameter. Note, these MRI 
parameters are based on enhancement kinetics and are independent of visual 
"morphology" cues used by the radiologist. Certainly, morphologic signs of 
smooth/rounded versus irregular/spiculated lesion borders are diagnostically important 
properties. Morphology-based reading of the MRI relative to mammography and the 
analysis, is more time consuming. This aspect of the study is ongoing using the existing 
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patient database. When completed, results of the clinical reading study will be 
summarized in a subsequent publication. 

RESULTS 

Spatial correlation between pathology, radiologic follow-up (>1 years) and ROIs 
quantified by MRI was achieved for 93 lesions/tissues with the quantitative results 
summarized in the Table below. 


Table 

Mean ± StdDev of MRI-derived Contrast Enhancement Kinetic Properties 


Tissue Type 

Enhancement 
Time Constant 
[sec] 

Enhancement 

Rate (~ P*S) 
riO -3 sec 4 ] 

Enhancement 
Amplitude 
[% Signal Change] 

Enhancement 
Amplitude x 
Enhancement Rate 

Benign (N=54) 

135 ± 364 

19 ±16 

132 ±139 

4 ± 11 

Malignant (N=39) 

27 ±40 

79 ±75 

211 ±101 

18 ±21 

p value 

0.052 

< 0.001 

0.03 

< 0.001 


Of the quantitative MRI parameters studies, enhancement rate provided the greatest 
discrimination between benign and malignant tissues in our database. Figure 2 illustrates 



0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 


FPF {1 - Specificity} 


ROC curves for Enhancement Rate, Amplitude and Rate x Amplitude. Good quality 
Xray mammography targets a Sensitivity >85% and Specificity >90%, and a positive 
predictive value (PPV) of 25%-40% for women recommended to biopsy [12]. The MRI 
ROC on quantitative (non-morphologic) data suggests a Sensitivity -85% and Specificity 
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=75%. The enhancement rate decision threshold at this Sensitivity/Specificity point was 
Enhancement Rate = 22/sec = 1000/Tau, which corresponds to an enhancement time 
constant of 43seconds. The MRIPPV was 72% at this decision threshold. The aggregate 
of individual tissue/lesion enhancement data are shown in Figure 3. 

Figure 3: Scatterplots of enhancement properties as a function of tissue type. Note, Tissue 
Types 1 - 5 are considered benign and Tissue Types 6 - 11 are malignant. (A) enhancement 
time constant; (B) enhancement amplitude; (C) enhancement rate, which is primarily 
dependent on PS product and (D) enhancement amplitude x rate. 


Figure 3A 
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Figure 3D 


Enhancement Rate x Amplitude 
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CONCLUSIONS 

Technical objectives of high temporal sampling and good spatial resolution for 
quantitative description of breast lesion enhancement properties were satisfied by the 
methodology developed for this study. Viewsharing increased temporal resolution to 5- 
6second sampling, however, a comparison of fitted kinetic parameters on viewshared data 
versus without viewsharing (i.e. 10-12second temporal sampling) suggest there is not a 
quantitative advantage to viewsharing. There is, however, a minor cosmetic 
improvement on cinegraphic images. Given this, we believe temporal sampling below 
lOsec imageset is not required for kinetic characterization of breast tissues. Moreover, 
newer MRI systems allow faster scanning thus 10-20second temporal sampling can be 
achieved without a substantial loss in spatial resolution, thus would be prefered over 
keyhole methods assuming the MRI system is capable to handle the data volume/image 
reconstruction. To date, we still prefer to perform these functions offline due to data 
handling restricions on our MRI systems. In addition, we believe motion correction prior 
to image reconstruction is desirable to reduced motion artifact in subtracted images. This 
too requires substantial post-processing overhead. Analysis of kinetic properties 
demonstrated that the rate of enhancement provided greates discrimination of malignant 
and benign tissues, although overlap of benign / malignant groups remain. These 
findings are consitent with those of other groups [11]. The amount (amplitude) of 
enhancement had very little diagnostic value. 

Kinetic properties can be combined with tissue/lesion morphology to further enhance 
specificity of breast MRI. Figure 4 illustrates an example where enhancement rate was 
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low (i.e. suggesting benign), but morphology shows an irregular lesion border which is 
suspicios for malignancy - this lesion was invasive ductal carcinoma by pathology. Thus, 
while morphology was not an emphasis in this project, clearly morphologic and kinetic 
properties are complimentary. 


Figtrre 4 




Figure 4: Irregular borders of this lesion suggest a malignancy even though its enhancement rate was low. 


In conclusion, high sensitivity of dynamic contrast-enhanced breast MRI, and the ability 
to detect additional lesions, even in dense breasts, suggests high spatial and high temporal 
resolution MRI methods will have genuine clinical value. Possible applications include: 
(1) survey breasts for additional foci prior to surgical biopsy of lesions detected by 
conventional means; (2) problem-solving imaging in dense tissues or those obscured by 
implant; (3) means of increased surveillance of women who "test positive" by future 
serologic or genetic breast cancer screens. 
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Dynamic Three-dimensional Imaging 
with Partial K-Space Sampling: Initial 
Application for Gadolinium-enhanced Rate 
Characterization of Breast Lesions 1 


PURPOSE: To evaluate a method to 
monitor gadolinium enhancement pat¬ 
terns at magnetic resonance (MR) imag¬ 
ing with high temporal resolution and 
full coverage through both breasts. 

MA TER1ALS AND METHODS: In 12 
patients with 13 masses, including 
nine carcinoma, nonenhanced three- 
dimensional MR imaging was per¬ 
formed with full-matrix resolution. 

At dynamic imaging, 32 serial passes 
were made during bolus administra¬ 
tion of contrast material, and tempo¬ 
ral resolution was reduced to 12 sec¬ 
onds by collecting the central (low 
spatial frequency) 32 x 16 or 16 x 16 
phase-encode views. Full-matrix dy¬ 
namic images were reconstructed by 
complementing central phase-encode 
data with precontrast data from pe¬ 
ripheral high-spatial-frequency views. 

RESULTS: Results at time-course 
analysis with a mono-exponential 
saturation model indicated malignant 
lesions tend to show rapid (<60 sec¬ 
onds) contrast change relative to be¬ 
nign masses and normal tissues. One 
cancer displayed an exceptionally 
slow contrast change (260 seconds). 

CONCLUSION: The technical objec¬ 
tives of full tissue coverage, rapid 
temporal sampling, and quantifica¬ 
tion of enhancement curves are met 
with this method for certain lesions 
(>5 mm in largest diameter). 

Index terms: Breast, MR, 00.12143 • Breast 
neoplasms, MR, 00.12143 • Magnetic resonance 
(MR), k-space 

Radiology 1995; 196:135-142 


E arly magnetic resonance (MR) 
imaging demonstrated only lim¬ 
ited success in assisting diagnosis of 
breast disease (1-4). Initially, only tis¬ 
sue-inherent differences in nuclear 
MR relaxation times, T1 and T2, were 
the primary sources of contrast. Al¬ 
though these contrasts are relatively 
high, it was observed that they have 
little value in helping detection and 
diagnosis of breast cancer (1,3,5,6). 
Variable admixtures of intense fat sig¬ 
nal and inadequate fat suppression 
have also seriously hindered image 
interpretation and detection of breast 
carcinoma. 

More recently, gadolinium-en¬ 
hanced MR mammography has gen¬ 
erated considerable interest (7-10). It 
is generally accepted that gadolinium 
enhanced MR imaging offers high 
sensitivity for the detection of breast 
cancers. The specificity of gadolinium- 
enhanced MR mammography to dis¬ 
tinguish benign from malignant tis¬ 
sues, however, has been limited by 
the fact that benign tissues can also 
exhibit contrast enhancement (9,11,12). 

The rate of gadolinium enhance¬ 
ment has also been studied as a po¬ 
tential discriminator of benign and 
malignant breast neoplasms (13-15). 

In their initial summary of findings in 
25 patients who underwent dynamic 
gadolinium-enhanced MR imaging, 
Kaiser and Zeitler reported that all 
malignancies displayed rapid en¬ 
hancement (ie, -100% increase in 
signal intensity within the first 2 min¬ 
utes after injection of 0.1 mmol/kg 
gadopentetate dimeglumine), whereas 
benign lesions demonstrated less and 
slower signal intensity change (13). 


These authors have since updated 
their dynamic gadolinium-enhanced 
MR imaging study, in which sensitiv¬ 
ity and specificity for cancer detection 
was greater than 95% (16). In other 
dynamic gadolinium-enhanced MR 
imaging studies, higher false-positive 
rates were observed; thus, there is not 
a consensus among investigators on 
methods or diagnostic value of this 
approach (11,12,17). Findings in these 
and other studies illustrate that fibro¬ 
adenomas and proliferative fibrocystic 
change, in particular, are sources of 
false-positive findings. Breast malig¬ 
nancies typically exhibit rapid en¬ 
hancement while benign lesions usu¬ 
ally, although not necessarily, exhibit 
slower contrast change. Malignant 
lesions may enhance slowly, thus re¬ 
ducing sensitivity (12,18). 

Use of various acquisition and en¬ 
hancement rate quantification schemes 
may contribute to discrepancies in 
findings (13-15). In assessment of 
enhancement rate, important param¬ 
eters include temporal resolution, 
volume of tissue covered, and the 
method used to quantify signal inten¬ 
sity change. In particular, important 
information on the rate of enhance¬ 
ment may be obscured if the temporal 
resolution is inadequate to sample 
rapid contrast changes. In addition, 
problems will obviously arise if the 
dynamic imaging volume includes 
an enhancing lesion only partially 
or does not include all lesions in the 
breast. These issues arise in part be¬ 
cause there is an inevitable trade-off 
in standard MR imaging between im¬ 
aged volume and temporal resolution. 

MR imaging with methods such as 
high spatial resolution, three-dimen¬ 
sional (3D) volumes, suppression of 
fat signal, and magnetization transfer 
contrast have been studied in an at- 


Abbreviations: ROI = region of interest 3D = 
three-dimensional 
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The time course of signal change in an 
ROl of' enhancing tissues was fit to a 
mono-exponential saturation model with 
use of a Levenberg-Marquardt least- 
squares optimization algorithm (23). Fitted 
parameters include contrast onset, t„; am¬ 
plitude of signal change. A; and enhance¬ 
ment time constant, t, used within the fol¬ 
lowing model: 


AS(/) = 


A[\ - e" ( '-' 0)/T ], for t > f 0 
0, for t < t {) . 


0 ) 


RESULTS 

Regarding mass conspicuity, con¬ 
ventional precontrast spin-echo Tl- 
and T2-weighted images frequently 
did not delineate location or spatial 
extent of the masses. Visual compari¬ 
son of pre- versus postcontrast images 
(Tl-weighted spin-echo and gradient- 
recalled-echo images) was adequate 
to identify nearly all (12 of 13) pal¬ 
pable and/or mammographically de¬ 
tected masses. Mass detection and 
localization were facilitated by inspec¬ 
tion of anatomic and subtracted cine- 
graphic images generated by means 
of the keyhole scheme. Subtraction 
removed many image features that 
partially obscured lesions, and cine- 
graphic display further delineated 
lesions in terms of time and intensity 
of contrast change. As the cinegraphic 
format provided the greatest conspi- 
cuity of enhancing tissues, it was used 
to interactively guide ROI definition 
on the time-series images. One super¬ 
ficial palpable mass (»1 cm in largest 
diameter) displayed signal intensity 
enhancement, although the enhance¬ 
ment was not distinct relative to sur¬ 
rounding tissues. The mass region 
was identified as a capsule on the skin 
that was visible at MR imaging. This 
mass was occult at conventional x-ray 
mammography and US and was 
shown to be benign at pathologic 
analysis. 

Figure 1 illustrates data obtained in 
a 56-year-old patient with pathologi¬ 
cally proved invasive lobular carci¬ 
noma in two adjacent sites separated 
by approximately 1 cm. Mammo¬ 
grams depicted only the larger mass. 
An early-time-point (ie, before con¬ 
trast intensity change) keyhole image 
(Fig la) does not depict the two ma¬ 
lignant foci, which were most appar¬ 
ent on the late-time-point subtraction 
keyhole image (Fig lb). Keyhole ma¬ 
trix size applied in this case was 256 x 
32 x 16. Several keyhole subtraction 
images from the 32-image time series 
for this section are shown in Figure 
lc. While the relatively rapid en¬ 
hancement period is visually appar¬ 
ent in the time series, temporal con¬ 


trast patterns of various tissues are 
best viewed graphically, as shown in 
Figure Id. Mono-exponential param¬ 
eters and fitted curves are also shown 
in this plot. The time constants of the 
lesions (t = 52 or 55 seconds) indicate 
that the majority of signal intensity 
change occurs in less than 1 minute. 
The number or location of these 
masses was not readily apparent on 
precontrast images (Fig la), although 
their region roughly correlated with a 
suspect mass seen on the mammo¬ 
gram. 

An example of a slowly enhancing 
lesion is Figure 2, which was obtained 
in a 56-year-old patient with a prob¬ 
able fibroadenoma. Findings at MR 
imaging and mammography were 
easily correlated for this particular 
mass, which had been sampled cyto- 
logically 7 years earlier. No change in 
size was seen, as noted in Materials 
and Methods. A subset of the subtrac¬ 
tion-image time series illustrates 
gradual enhancement of the lesion 
(Fig 2c). The keyhole matrix size ap¬ 
plied in this case was also 256 x 32 x 
16. As with the data in Figure Id, the 
time-course curve is well fit by an ex¬ 
ponential (Fig 2d), although for this 
lesion the time constant is substan¬ 
tially longer (t = 119 seconds). 

A summary of fitted mono-expo¬ 
nential time constants for a variety of 
tissue types and all patients is shown 
in the scatter plot of Figure 3. ROIs 
from tissues labeled "normal" were 
derived from areas remote to the mass 
that displayed sufficient enhancement 
for curve fitting but were otherwise 
normal in appearance. These were 
typically regions of diffuse enhance¬ 
ment from either the ipsilateral or 
contralateral breast. No more than 
one normal-tissue ROI was derived in 
each patient. Two of the five invasive 
ductal carcinoma tumors were found 
in one patient (within one breast and 
separated by approximately 5 cm), 
and two of the three invasive lobular 
carcinoma tumors were found in one 
patient (Fig 1). ROIs of tissues adja¬ 
cent to invasive ductal carcinoma tu¬ 
mors were defined to avoid signal 
blur from the lesion focus. As men¬ 
tioned in Materials and Methods, the 
keyhole matrix size was 256 x 16 x 16 
in three of 12 patients. Masses charac¬ 
terized in these patients were three 
invasive ductal carcinoma tumors 
(with largest diameters of 18,13, and 
22 mm and t of 14, 20, and 18 seconds, 
respectively) and one neuroendocrine 
carcinoma (largest diameter, 14 mm; 
t, 37 seconds). 

There was a clear tendency for the 
malignant lesions to display rapid 


contrast changes (t < 60 seconds) 
relative to benign or normal tissues. 
However, one invasive lobular cancer 
displayed an exceptionally long time 
constant with t of 260 seconds. Patho¬ 
logic examination revealed the tumor 
had a largest diameter of at least 1 cm 
(margin was positive at initial biopsy), 
gross size of 1.2 x 1.0 x 0.8 cm, no 
angiolymphatic invasion, and no ex¬ 
tensive ductal carcinoma in situ; in 
addition, 17 of 17 lymph nodes were 
negative for malignancy. 

DISCUSSION 

We tested a technique with which 
to perform dynamic, contrast materi¬ 
al-enhanced MR imaging of the en¬ 
tire volume of both breasts with a 
temporal resolution of approximately 
12 seconds. The usefulness of the 
method was explored in our study of 
12 patients with a total of 13 breast 
masses. As found in other studies, we 
noted a trend for malignant cancers to 
display rapid enhancement relative to 
benign and normal tissues, although 
the presence of one slowly enhancing 
cancer was observed. This would be 
considered a false-negative case if the 
enhancement time constant were 
used as the sole diagnostic indicator 
of malignancy. The small number of 
patients in this study, however, pre¬ 
vents us from drawing statistically 
meaningful inferences regarding the 
diagnostic value of this technique. 
Technical issues aside, other investi¬ 
gators have noted both false-negative 
and false-positive findings relative to 
enhancement rate (11,12,17,18). 

The viability of this approach for 
management of breast disease neces¬ 
sitates determination of sensitivity 
and specificity rates for each applica¬ 
tion. The role that we are currently 
evaluating is the contribution of in¬ 
cremental information beyond that 
provided with conventional imaging 
modalities when a suspect mass is 
detected. In this capacity, emphasis is 
placed on characterization of previ¬ 
ously detected masses that are usually 
more than 5 mm in largest diameter. 
Lesion localization for eventual dy¬ 
namic gadolinium-enhanced MR im¬ 
aging is hindered by nonspecific le¬ 
sion properties in precontrast MR 
images and/or by poor correlation of 
findings with those of other imaging 
modalities, owing to variable breast 
geometries and views. This problem is 
essentially solved with the 3D key¬ 
hole approach, which has the inher¬ 
ent features of complete coverage 
through both breasts and high tempo¬ 
ral sampling of dynamic contrast 
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merit properties. An alternative to the 
keyhole method is rapid acquisition 
of low-spatial-resolution data sets that 
are reconstructed and viewed inde¬ 
pendently from high-resolution data 
sets. However, we believe the com¬ 
plexities of keyhole reconstruction are 
justified to ensure the temporal en¬ 
hancement profiles and spatial prop¬ 
erties of a region are correctly related. 
Spatial blur of enhancement that re¬ 
sults from acquisition of limited k 
space will spread along phase-encode 
directions to bordering tissue. That is, 
truly unenhancing tissues can dem¬ 
onstrate apparent enhancement that 
is actually caused bv bleeding of sig¬ 
nal intensity from adjacent enhancing 
tissues. Vertical smearing of signal 
intensity from small enhancing ob¬ 
jects is illustrated in the subtraction 
images of Figures lb and 2b. 

In addition, the apparent amplitude 
• *f enhancement within trulv en¬ 


A Figure 2. Data from a coronal image obtained through the breast of 
~ 56-year-old patient with a stable fibroadenoma, (a) Early-time-point (ie, 
precontrast) keyhole image, (b) Subtraction image from a late time point 
(t = 355 seconds) of dynamic imaging shows contrast enhancement in 
the mass (ROD. (c) Subset of time series for this section illustrates a grad¬ 
ual contrast change. " = seconds, (d) Time-course analysis of the mass 
indicates that this fibroadenoma has a moderate to long time constant 
(t = 119 seconds). 


The apparent enhancement time 
course of the lesion is derived by mul¬ 
tiplying by an ROI mask and 

then performing spatial integration. 
The convolution theorem allows us to 
relate this in k space as 

AROI(f) = ROI(f) - ROI(O) 

-cw/ikM 5 *) 

• L(k y ,k z ) 

• K(k y ,k z )dk y dk z/ (4) 

where R(k y ,k z ) is the Fourier transform 
of the ROI mask. If we assume the spa¬ 
tial pattern of the ROI closely matches 
that of the enhancing lesion, then 
R(ky,k z ) = L(kyk z ); thus, the integrand 
becomes | L( y ,k z ) p within the keyhole 
area and zero elsewhere. Finally, the full 
matrix acquisition is used to normalize 
the contrast suppression effect: 

AROI(f) = P -C(f) • ROI(0); 

where 


hanced tissue will be diminished. This 
effect in the phase-encode dimen¬ 
sions, y and z, is demonstrated with 
the function that characterizes a uni¬ 
formly enhancing lesion as 

L(y,z,t) = L(y,z,0) • [1 + C(t)l 
= .7/- , |L(k y ,M) 

• [1 + C(t)]\, (2) 

where CO) defines the temporal prop¬ 
erties with C(0) of 0 and L(k y ,k z ,0) de¬ 
scribes spatial properties of lesions in 
k space. Changes in contrast in the 
reconstructed keyhole image of this 
lesion, _/(y,z,t), are limited by the 
spatial frequencies that lie within the 
sampled keyhole matrix. Thus for 
K(k v k,) = 1 within the keyhole range 
and zero otherwise, we have 

y (y,z,t) = / 7 -l |L(k y ,k z ,0) 

(1 + KfMO • C(f)]|- (3) 


0 “ .keyhole area IUMO1 2 dkydk* 

+ JfuU k space I L(ky,k z ) | 2 dkydk z | (5) 

The factor p represents the degree 
of suppression in contrast change that 
results from sampling only a portion 
of k space. For the above derivation to 
apply, the time to complete each key¬ 
hole pass must be short or must be at 
least comparable to the period of 
strong contrast change. The lesion 
contrast state is considered nearly 
constant for each keyhole pass but is 
allowed to vary from pass to pass. 
While coverage of full k space may 















Figure 4. (a) Maximum intensity projection of subtraction images in the superior-inferior projection. Six noncontiguous time frames are 

shown from images normally viewed cinegraphically. ROm delineates an invasive ductal carcinoma. Enhancement in linear structures delin¬ 
eates vessels. A = anterior. L = left. P = posterior, R = right, f = time (seconds). " = seconds, (b) Correspond, ng enhancement-hme^course 
curves for ROIs 1-f illustrate changes in relative rate and intensity of enhancement. Curves were derived from ROIs drawn on related single- 
section cinegraphic displays. ROl# I = left ventricular blood (broken line), RO/#2 = liver (O), RO/#3 - invasive ductal carcinoma (solid line), 
ROI#4 = area adjacent to mass (dashed line). 


across these studies. The immediate 
technical objectives of this study were 
to overcome these critical limitations. 
These objectives are largely met with 
the 3D keyhole technique, which of¬ 
fers full tissue coverage at high tem¬ 
poral sampling for subsequent quantita¬ 
tive enhancement time characterization. 
As we have discovered, the ability to 
dynamically study entire breasts at MR 
imaging, thus obviating prospective 
identification of one or all lesions amidst 
complex contrasts, is a great practical 
advantage. A drawback of this approach 
is the large data volume generated per 
case. Currently, we resolve the data 
management issue by performing recon¬ 
struction, display, and analysis on an 
off-line workstation with customized 
software tools. Our software has been 
specifically designed to present data in a 
variety of formats for specific purposes, 
such as the survey of all tissues in cine- 
graphic maximum-intensity-projection 
images of subtraction images or the cine- 
graphic viewing of a single section for 
quantitative analysis of individual tis¬ 
sues. This method also benefits from the 
other well-known features of breast MR 
imaging, such as the ability to image 
breasts with dense parenchyma or im¬ 
plants and to provide high-quality refor¬ 
matted tomographic views that can po¬ 
tentially aid planning for surgery. ■ 
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Quantitative analysis of dynamic gadolinium-DTPA^enhanced MR|is emerging as a, h ighly wici tiva ydifoA' * 


tool for detecting malignant breast tissue. Three-dimensiond"fapT(runagilfg^e6hniaues, such as 
keyhole MRIJ^BM^petj£^«§gMQp«^p*aglsgV high temporal sampling rates to accurateljmack 
contrast enhancement and washout in lesions over the course of multiple volume acquisitions?' 
Patient motion during the dynamic acquisitions is a limiting factor that degrades the image quality, 
particularly of subsequent subtraction images used to identify and quantitatively evaluate regions 
suggestive of malignancy. Keyhole imaging is particularly sensitive to motion since datasets ac¬ 
quired over an extended period are combined in k-space. In this study, motion is modeled as set of 
translations in each of the three orthogonal dimensions. The specific objective of the study is to 
develop and implement an algorithm^) correct the consequent phase shifts in k-space data prior to 
offline keyhole reconstruction of imagfc of M> volume breast MR acquisitions. © 1999 American 
Association of Physicists in Medicine. [S0094-2405(99)01005-6] . . _ 1 
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INTRODUCTION 

Contrast enhanced MRI is a developing technique for char¬ 
acterization of breast lesions. 1 " 8 This method is based on the 
tendency of lesions to exhibit enhancement after administra¬ 
tion of gadolinium-based contrast media. Furthermore, ma¬ 
lignant tissues typically exhibit rapid enhancement due to 
anomalous vascularity, in which the combination of in¬ 
creased capillary density and permeability results in rapid 
enhancement, and occasionally washout, of the contrast 
agent in a dynamic contrast MRI study. 9 " 12 Benign tissues, 
on the other hand, tend to have a slower enhancement rate. 5 
Some benign tissues, such as fibroadenomas, have been ob¬ 
served to exhibit rapid enhancement, 13 thereby reducing the 
specificity of dynamic contrast enhanced MRI of the breast. 
Nevertheless, analysis of the temporal enhancement patterns 
is considered an emerging tool that can be used for detection 
and differentiation of various breast lesions. For this pur¬ 
pose, two-dimensional multisection breast MR imaging has 
been performed with an associated temporal resolution of 
~1 min. 5 Even better temporal resolution (15 s) has been 
achieved by imaging one or a few pre-selected slices. 7 Real¬ 
istically, selection of a few slices based on findings at x-ray 
mammography, ultrasonography, or physical examination 
can be difficult and is useful only in characterizing known 
lesions. Correlation of these findings to a few MRI imaging 
slices is difficult due to the variable deformation of the 
breast Also, since certain lesion are not easily identified in 
the absence of contrast agents, a precontrast MR image by 
itself could prove inadequate for prospective selection of ap¬ 
propriate imaging sections. Thus for dynamic MR breast im¬ 
aging to be a useful screening tool, techniques that provide 


3D volume images of one or both breasts while sampling 
contrast enhancement information at high temporal resolu¬ 
tion, arc desired. 

Three-dimensional dynamic MRI consists of acquiring se¬ 
rial 3D volume images over a finite time period that encom¬ 
passes contrast enhancement. The time per acquisition is a 
function of the imaging parameters, such as matrix size, 
number of excitations (NEX), and the type of acquisition 
sequence used. 14 Three-dimensional volume imaging with 
fat suppression using magnetization transfer contrast has 
been performed to yield images of good anatomic quality, 
i.e., high spatial resolution, but with acquisition times rang¬ 
ing to several minutes. 15 Rapidly enhancing tissues have 
been observed to reach peak enhancement on the order of 
tens of seconds. To be sensitive to the most rapidly enhanc¬ 
ing tissues high sampling rates are required which is 
achieved, in part, by using fast imaging techniques, such as a 
3D fast rf-spoiled gradient recalled echo (3D SPGR) se¬ 
quence used in this study. Despite improvements in the 
speed of the imaging pulse sequence, there remains a 
tradeoff between high spatial resolution and high temporal 
sampling rates, which relates to the objective of breast lesion 
detection and characterization. One approach to address this 
tradeoff has been to subsample A-space using a tailored tra¬ 
versal pattern, while maintaining a high dynamic temporal 
sampling rate. The resulting spatial resolution usually de¬ 
pends upon the extent of A-space traversal and the recon¬ 
struction algorithm. Examples include reduced-encoding im¬ 
aging using generalized reconstruction (RIGR), 16 use of 
projection reconstruction trajectories, 17 and dynamic spiral 
breast imaging. 18 A particular technique, used for data acqui¬ 
sition in this study, that similarly mitigates the spatial- 
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temporal tradeoff is the “keyhole” MRI technique described 

’ in* the following section. 

„ The keyhole technique involves the acquisition of one or 
more reference high-resolution datasets followed by a dy- 
namic series of datasets having a limited extent in Jk-space. 19 

• Tne reduced matrix acquisition usually samples only the cen¬ 
tral portion of k-space, which is then spliced with peripheral 
k-space data from a full matrix reference dataset acquired 
initially. In keyhole acquisitions it is assumed that the sig¬ 
nificant dynamic events (such as amplitude changes in breast 
lesions), are tracked by acquiring the central k-space lines. 
The adequacy of tracking changes depend on the spatial fea¬ 
tures of the enhancing object. It is an inherent limitation of 
keyhole imaging that features exhibiting temporal enhance¬ 
ment will be blurred by an amount dependent on the keyhole 
matrix size and the spatial dimensions of the object. Fine 
structures with a higher peripheral spatial frequency content 
will be blurred by a greater amount. 20 However, the spatial 
detail of static, nonenhancing, background structures is re¬ 
stored by splicing together the high spatial frequencies from 
the reference dataset. This serves to provide a high- 
resolution anatomic context for evaluating the dynamically 
changing lesions. 21 

A significant source of error related to keyhole imaging 
will be due to breast motion over the total acquisition period. 
This will manifest as phase differences between the reference 
and dynamic datasets. For contrast enhanced MRI, quantita¬ 
tive analysis is commonly performed on subtraction images 
that emphasize temporal changes. Keyhole reconstruction of 
the phase-deviant datasets will result in substantia] edge ar¬ 
tifacts and blurring in the subsequent subtraction images 
used for quantitative analysis. 22 A number of motion models 
and corresponding post-processing techniques have been 
proposed to reduce motion artifact for two-dimensional Fou¬ 
rier transform imaging. 23-27 The purpose of this paper is to 
outline a 3D model for motion during the rapid dynamic 
acquisition and describe a method of estimation and correc¬ 
tion for the phase artifacts introduced by motion. The motion 
correction algorithm was initially verified by computer simu¬ 
lation. It was tested on a phantom experiment dataset and 
applied to clinical breast studies. 

MATERIALS AND METHODS 
Motion correction 

Dynamic breast imaging data was acquired using a dual 
phased array coil, which permitted processing of data from 
the right and left coils separately. Thus motion correction 
was applied individually to each breast. The motion model 
was based on a consideration of the patient configuration and 
imaging rates. With some mild compression most patients 
were reasonably well constrained within the breast coil. The 
3D SPGR imaging sequence acquired a single 3D volume 
once every ~10 s. Since the patients were cautioned to hold 
still, the most likely causes of motion were gradual posi¬ 
tional shifts that evolved on a time scale longer than the 
dynamic temporal sampling rate. Thus, we assume negligible 
intrakeyhole motion and consider a mean position over each 


dynamic volume acquisition. Based on visual inspection of 
the subtraction edge artifacts observed, we further limit our 
analysis to simple translations in each of the three orthogonal 
dimensions. Our primary focus, therefore, was to detect and 
correct simple displacements that occurred over the duration 
of the overall dynamic acquisition on a per keyhole basis. 
Thus, the model assumed that after the reference acquisition 
and between each subsequent dynamic acquisition, die indi¬ 
vidual breasts were allowed to move independendy as rigid 
bodies undergoing only transladonal motion. 

It is known that a simple displacement in space introduces 
a corresponding phase shift in the spatial frequency signal 
while the magnitude of the data remains unchanged. 28 Thus 
there is a phase difference between the reference dataset and 
the dynamic dataset acquired after the object has been dis¬ 
placed. Let S(k x ,k y ,kj) represent the reference dataset in the 
spatial frequency domain, given by 

S(k x ,k yt k t )= f p(x,y,z)‘ Wk ‘ x+k >> +k S»dxdydz. ( 1 ) ' 

Now, let the object, p(x,y,z), be displaced by Ax, Ay, Az 
at a given Ume point. The new posidon of the object can be 
described as a convolution with a displaced delta function 

p'{x,y,z)=p(x,y,z)® S(x- Ax.y - Ay ,z - Az). (2) 

Consequently, the k-space data from the object is now given 
by 

S'(k x ,k y ,k z )=S(k x ,k y tkt yl-‘«A*+V>+ k zW). (3) 

Thus the net phase shift due to 3D translation is given as 

<(>{k x ,ky,k t )=k x &x+kyby+k t hz. (4) 

In keeping with our motion model, the objective of the 
motion correction algorithm was to estimate this linear phase 
shift of the entire keyhole k-space block relative to the ref¬ 
erence data block per spatial frequency axis, and apply a 
phase correction to each dataset prior to keyhole reconstruc¬ 
tion. The approach was to compute an averaged linear phase 
roll along each individual k-space axis, [Eq. (5)] which is in 
effect a projection of mean phase difference over the 3D 
dataset, onto that particular axis. 

To implement the algorithm, first a phase difference ma¬ 
trix was generated between a central (32X32X32) kernel 
extracted from the reference dataset and each dynamic time 
point The phase difference was estimated on reduced matrix 
spatial frequency datasets, because it cannot be reliably com¬ 
puted in the low amplitude, peripheral k-space, regions of the 
signal. In these regions the noise amplitude is comparable to 
that of the signal and hence the phase could fluctuate ran¬ 
domly, obscuring the effect of linear motion. The average 
phase difference along each spatial frequency axis was com¬ 
puted as follows: 

*<**>= 2 kyl, k-d>(k x ,ky,k t ) 

(n y ‘n t ) -14/ 

+constant term, 
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m r 'I/ 

’ 4 + constant term, (5) 

^ k *t~7n -it \ k y'M k x> k y’ k i) 

K x y) ^ l 

+constant term, 

where, n X9 n y9 n z , is the matrix size and, <f>(k x ) 9 <f>(k y ) 9 
<f>(k y ), are the averaged phase deviations per &-space axis. 

Thus we obtain a phase difference function per &-space 
axis which istlneasure of the averaged phase difference be¬ 
tween the reference and dynamic datasets, due to displace¬ 
ment in the corresponding spatial axis during a particular 
dynamic timepoint. For example, in Fig. 1(a) the coronal 
plane image shows a noticeable edge artifact. Figure 1(b) in 
turn, shows a net averaged linear phase roll detected along 
the k x and k z dimensions and negligible phase deviations in 

k y . This implies linear displacement along the x and z di¬ 
mensions. Thus the phase function is an estimate of the 
translation induced phase shift. The linear phase function 
could contained wrapped phase depending on the extent of 
motion along a particular spatial axis. Thus, an unwrapping 
algorithm was applied to each function so that it increased or 
decreased monotonically. These averaged phase deviation 
functions were submitted to a linear least squares estimation 
algorithm from which a representative slope per spatial fre¬ 
quency axis was computed. This slope was used to correct 
the phase deviant dataset by simply applying an inverse lin¬ 
ear phase ramp such that 

^cotr(^ 4 c yky,k z ) k x s x "h k y s y ^~k z s z , (6) 

where, s x9 s y9 and s z were the fitted slopes in k x9 k y9 and 

k z . 

The constant phase shift term between the reference and 
the slope-corrected datasets was also determined and incor¬ 
porated in the phase correction algorithm. In the case of very 
large displacements [>£ (FOV) field of view] along each 
axis, the phase difference function could be wrapped twice. 
The unwrapping technique used would not be able to account 
^^j^this in a single stage and hence the correction algorithm 
was applied to the data in two consecutive stages. 

ANALYSIS OF AMPLITUDE MODULATION 
EFFECTS 

The amplitude modulation of the data that is associated 
with uptake of Gd-DTPA could theoretically, be misinter¬ 
preted as a phase shift Computer simulations were used to 
further analyze this effect We investigated the effect of sig¬ 
nal enhancement by considering signal increases starting 
from 100 to a maximum of 500 percent Two simplified ana¬ 
tomic configurations of breast fat and lesion were consid¬ 
ered. First we varied the percentage of simulated lesion area 
to breast fat area from l%-60%. We increased lesion size 
relative to the size of the background fat while changing the 





Fig. 1. (a) Reconstructed image of single (right) breast after separating dual 
phased array data. Bottom image shows representative subtraction artifact 
on coronal slice of dynamic timepoint 14. (b) Averaged phase deviations per 
spatial frequency axis for dynamic timepoint shown in (a), relative to the 
reference dataset. 
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percent signal enhancement within the lesion. The phase dif¬ 
ference between each enhanced and the original unenhanccd 
dataset was computed. Second, we varied the displacement 
between the axes of the simulated breast fat and lesion from 
0 to 5 cm. We displaced breast fat relative to a fixed lesion. 
Again, for each axial displacement configuration the percent 
signal enhancement was varied and the phase difference be¬ 
tween the enhanced and unenhanced datasets was computed. 

MR imaging 

The breast studies were performed on a 1.57 system 
(General Electric Medical Systems, Milwaukee, Wisconsin), 
using a dedicated breast dual phased array coil. For the clini¬ 
cal data, 71 and 72 weighted axial scans were performed 
initially. A pre-contrast coronal reference dataset was ac¬ 
quired using the locally developed 3D fast rf-spoiled 6KB- 
(3D SPGR) sequence, 40 degree flip angle. Imaging param¬ 
eters were a TR/TE of —13/5 ms, acquisition matrix, 256 
X128X32 with 3-5 mm thick sections and four excitations. 
The field of view (FOV) ranged from 28 to 36 cm (typically 
32 cm). The dynamic contrast-enhanced series consisted of 
20 serial 3D volumes acquisitions with a reduced matrix of 
256X32X32 and single excitation, followed by one full ma¬ 
trix dataset. The dynamic segment spanned —5 min during 
which a bolus injection of Gd-DTPA was administered 
within the first 30-45 s of the scan, at a dose of 0.1 mmol 
per kilogram body weight. Each dynamic 3D dataset was 
acquired at 12.3 s intervals, which is equivalent to the overall 
temporal resolution. Finally a series of post-contrast coronal 
3D SPGR and an axial 71 weighted scan were performed. 

Keyhole reconstruction of the data was performed offline 
using custom software programs developed in AVS (Ad¬ 
vanced Visual Systems Inc., Waltham MA). The low spatial 
frequencies from each dynamic dataset were spliced into the 
reference dataset to create the corresponding full matrix dy¬ 
namic dataset, for subsequent 3D Fourier transform recon¬ 
struction. Cinegraphic loops of reconstructed anatomic and 
subtraction images were reviewed on the workstation. A pre¬ 
contrast time point was used as the subtraction mask. 



24.*" 73.*" 123" I72J" 221.4" 

Fig. 2. Coronal slice of 3D volume acquisition (a). Anatomic image of 
uncorrected coronal slice, (b) Subtraction image of the same slice with 
clearly visible subtraction edge artifact, (c) Time series of edge artifact in 
the same slice prior to motion correction, 49.2 s intervals, (d) Time series of 
edge artifact in the same slice after motion correction, 49.2 s intervals. 

displacement commonly observed over the duration of the 
scan. Uptake of contrast in a lesion was simulated by injec¬ 
tion of 10 cc of 2 mmol Gd-DTPA into 5 cc vials located in 
each breast phantom compartment 

Patient studies 

The clinical data included 64 breast patients who were 
scanned as part of an ongoing research project Study popu¬ 
lation consisted of patients with abnormal mammograms or 
ultrasonically detected suspicious masses who were referred 
for MR scans. Patients were scanned prone with their breast 
suspended in the coil. Mild compression was applied to the 
breasts but no rigid immobilization techniques were em¬ 
ployed. The standard clinical procedure was to caution the 
patients to hold still during and in between scans. 


Phantom studies 

The motion-correction algorithm was calibrated experi¬ 
mentally. Displacements were introduced in each of the three 
spatial axes using a lever arm rigidly attached to the experi¬ 
mental phantom located in the breast coil. Motion was in¬ 
duced by translating the lever arm over a calibration scale 
such that the phantom was displaced by an exact amount, in 
increments of 2 mm. The maximum displacements intro¬ 
duced were 2 cm along the x and y axes, and 1 cm along the 
z axis. 

The 3D dynamic simulation experiment was performed 
on a breastTinimicking phantom, on the 1.57 GE system. 
Mineral oil was used to simulate breast fat and water to 
simulate breast parenchyma in each of the phantom compart¬ 
ments. Motion was induced in only one breast phantom com¬ 
partment by raising and laterally displacing a lever attached 
to the phantom. Motions were designed to simulate the slow 


Quantitative assessment 

The effectiveness of the motion correction algorithm was 
quantitatively assessed by comparing subtraction edge arti¬ 
fact prior to and after motion correction. A single measure 
that summarized edge artifact was devised and computed by 
the following procedure. For each of the 64 breast cases a 
representative mid-breast slice wa s chosen. Based on the 
subtraction images, a single ROrencompassing all of the 
noticeable motion artifact induced edge was defined (Fig. 2). 
This ROI was used as a mask over the entire sequence of 
dynamic timepoints and a mean ROI intensity per time point 
was computed, for the corrected and uncorrected datasets at 
the selected slice of interest. Next a baseline value equal to 
the mean ROI intensity from the second time point in a series 
was subtracted off from all other time points, providing the 
desired estimate of subtraction edge artifact. This temporal 
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series of mean artifact ROI’s was further approximated to a 
■ <$inglo mean and peak artifact ROI over the entire dynamic 
acquisition. Similar mean ROFs intensities were calculated 
. f° r a region in each dynamic time point that corresponded to 
thd background noise signal, for both corrected and uncor- 
rec w ted datasets. In order to compare the edge artifact correc¬ 
tion across patients, the artifact ROI intensities were normal¬ 
ized by the mean background noise ROI intensity. 

RESULTS 

The calibration experiment for the motion correction al¬ 
gorithm yielded the following results. The correlation coef¬ 
ficient between the induced and detected displacements was 
0.999 for the x t y, and z dimensions. Thus, there was a close 
overall agreement between the induced and detected dis¬ 
placement along each of the spatial axes. A minimum dis¬ 
placement of 2 mm was induced and detected along each 
spatial axis. The theoretical limit to the maximum displace¬ 
ment that could be corrected for, is equivalent to a half field 
of view along that particular spatial axis. 

The phase slopes in k xt k y , and k z measured in the 3D 
dynamic experimental data, corresponded to a maximum dis¬ 
placement of —6 mm along the anterior-posterior direction 
(A/P), 2 mm along the right-left (R/L) direction and negli¬ 
gible motion along the superior-inferior (S/I) direction. Fig¬ 
ure 3 demonstrates the effect of motion correction on the 
phantom. A marked reduction in motion-induced blurring 
and ghosting was observed on the subtracted and nonsub- 
tracted images. An overall improvement in the edge artifact 
at breast boundaries as well as good background suppression 
was achieved for the motion-corrected subtraction images. It 
is interesting to note that after motion correction the vial on 
the right [Fig. 3(d)] is no longer visible. This is because the 
vial is truly outside the level of the selected slice in the 3D 
volume. It is visible on the uncorrected image due to sub¬ 
traction artifact. After motion correction the slice is reregis¬ 
tered to the appropriate cross sectional level. 

There was no significant difference in the mean ROI com¬ 
puted for each vial due to motion correction compared to the 
dataset without motion correction. This result was expected 
since the size of the simulated lesion was fairly large and 
therefore keyhole imaging could accurately track dynamic 
changes. 

The simulation investigating the effect of contrast modu¬ 
lation on phase estimation showed enhancement effects to be 
negligible. We found that, as long as the high signal unen¬ 
hancing breast fat, was larger in area or mostly coaxial with 
the enhancing lesion, the phase deviations were very small 
for all conditions of enhancement within the lesion. The 
computed phase errors were negligible for lesion as large as 
60% of breast area with axial separations up to 5 cm. This 
was validated by the phantom experiment where the lesion- 
simulating vial was small compared to the breast fat com¬ 
partment and located coaxially. There was no significant dif¬ 
ference between the mean ROI in the signal enhancing vial 
without motion, and the post-correction mean ROI in the 
signal enhancing vial with induced motion. This implies that 
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Fig. 3. Phantom experiment with breast mimicking phantom, coronal slice 
of 3D volume, (a) Unsubtracted original image, (b) Unsubtracted image 
after motion correction, (c) Subtraction image prior to motion correction, (d) 
Subtraction image after motion correction. 


the motion-correction phase estimates were computed accu¬ 
rately in the contrast-enhanced displaced datasets. Addition¬ 
ally, in a couple of patient datasets we observed substantial 
enhancement in the breast parenchyma with contrast admin¬ 
istration, however motion correction did not result in a mis- 
positioning of this enhancing parenchyma, since it is mostly 
coaxial with breast fat. 

There were a range of results for the clinical breast data, 
depending on the nature of the motion that might have oc¬ 
curred during the exams. For a typical study we maintained a 
32 cm FOV with 4 mm thick slices at a (256X32X32) ma¬ 
trix size. At these settings the maximum detected slopes, 
over all patients, corresponded to a displacement of 8, 5, and 
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• 3 mm in A/P, S/I, and R/L directions, respectively. 

* Qualitative visual assessment of the subtraction images 
r showed that in certain cases there was substantial observable 
* ^ motion artifact in the uncorrected datasets, that was success¬ 

fully removed in the corresponding motion-corrected images. 
v The improvement was most noticeable as a marked reduction 
in subtraction edge artifact. Low amplitude ghost artifact and 
blurring in the background was also substantially reduced 
(Fig. 4). This effect was improved clarity of 

features such as the nipples. Motiorrcorrection also enhanced 
the visual quality of fine structures such as vessels seen in 
axial and cross-sectional orientations (Fig. 2). From the plots 
of mean slope per time point in Fig. 4, we see that the motion 
correction algorithm has the freedom to individually correct 
each breast. In this example, greater displacements were de¬ 
tected for the right as compared to the left breast, providing 
the same overall degree of correction for both breasts. 

Quantitative analysis of motion artifact reduction for the 
clinical cases is summarized in Fig. 5. The mean subtracted 
artifact ROI intensity over all dynamic timepoints, normal¬ 
ized to the background noise signal intensity, motion cor¬ 
rected vs uncorrected datasets, is plotted in Fig. 5(a). The 
corresponding plot for the peak artifact ROI is shown in Fig. 
5(b). The line of unity represents equivalent artifact in mo¬ 
tion corrected and noncorrected datasets. The filled circles in 
the graphs [Figs. 5(a) and 5(b)] correspond to the peak and 
mean artifact correction levels for the edge artifact seen in 
Fig. 4. For this representative artifact we see that there was a 
roughly 40% reduction in subtraction edge artifact, resulting 
in the improved visibility described earlier. We can summa¬ 
rize the results of motion correction in terms of the edge 
artifact measure as follows: 

Motion correction produced a reduction in mean artifact 
ROI in a number of cases, i.e., there was an improvement in 
edge artifact suppression. For certain cases there was no sig¬ 
nificant difference between the mean ROI’s for the corrected 
and uncorrected datasets. This suggests that other sources of 
phase artifact, such as rotations and distortions, that did not 
fit the linear three-dimensional translational model, could be 
present. In none of the cases was the mean artifact ROI 
greater for the motion corrected dataset compared to the un¬ 
corrected dataset This implies that the correction algorithm 
did not introduce any additional artifact. 

CONCLUSIONS AND DISCUSSION 

Contrast enhanced MRI is developing into a fairly useful 
tool for the detection and characterization of breast tumors. 
The potential for breast MRI as a clinical diagnostic tool lies 
in the ability to achieve volume imaging of both breasts at 
high spatial resolution, yielding good anatomic detail. In ad¬ 
dition contrast enhanced dynamic imaging provides func¬ 
tional information that could assist in tumor characterization. 
The fidelity of the MRI image data is often limited by artifact 
sources including motion during the acquisition. 

In keyhole substitution MRI an effect of linear motion is 
to introduce phase discontinuities and shifts between the dy¬ 
namic and reference datasets. These phase deviations mani- 
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Fig. 4. Coronal slice of 3D volume acquisition, (a) Subtraction image before 
motion correction, (b) Subtraction image after motion correction. Note im¬ 
proved clarity of nipples (arrow) and reduced subtraction edge, (c) Mean 
slopes in kx, ky, and kz. per time point in left breast, (d) Mean slopes in kx, 
ky, and kz, per time point in right breast 


fest as edge artifacts and cause blurring, rather than an ob¬ 
servable gross displacement, of small lesions in subtraction 
images that are used for quantitative analysis. This is be¬ 
cause the bulk (three fourths) of the spatial information is 
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Motion Corrected VS Non Motion Corrected: Mean 
* Subtracted ROI Normalized to Background Noise Level 




Fig. 5. (a) Comparison of mean artifact level per dynamic series between 
motion-corrected and uncorrected datasets, (b) Comparison of peak artifact 
level per dynamic series in motion-corrected and uncorrected datasets. 


derived from the peripheral high frequencies In the reference 
dataset Figure 6 is an illustration of this effect. We observe 
that motion produces subtraction edge artifact and some 
amount of lesion blurring, but retains most of the lesion con¬ 
trast information. The artifacts can, however, obscure lesions 
or result in reduced visibility in subtraction images that nor¬ 
mally offer the greatest lesion conspicuity. Superficial le¬ 
sions that lie near high contrast boundaries could remain 
undetected. Furthermore, identifying lesion ROI*s in these 
artifactually blurred images could lead to inaccurate quanti¬ 
fication of contrast uptake characteristics. Thus, effective 
motion correction algorithms are desirable in contrast- 
enhanced breast MRI. 

The scope of this study was limited to gross translational 
rigid body motion of the breast The motion correction algo¬ 
rithm assumes no intrakeyhole motion, i.e., motion during 
the acquisition of a single dynamic block. In this* context, 
calculating averaged phase slope tends to identify gross 
shifts between the reference and each subsequent dynamic 
acquisition. Further, averaging permits the algorithm to be 
less sensitive to random phase fluctuations that it would be 
highly susceptible to, if we were to perform a local phase 



(b) 


Fig. 6. Example of subtraction edge due to keyhole reconstruction of motion 
corrupted dataset, while contrast information is retained, (a) Unconected 
coronal slice, showing malignant lesion and broad subtraction edges, (b) 
Same slice after correction, showing reduction in edge artifact, while main¬ 
taining lesion clarity and contrast. The detected displacements were 1.6,4.7, 
and 8.1 mm along the ac, y, and z axes. 


correction per point in *-space. However, there are other de¬ 
grees of motion such as rotation, distortions and respiratory 
and cardiac motions, that are commonly encountered during 
the scan. The assumption of no intrakeyhole motion itself 
does not strictly hold and this combined with heart motion 
contributes to the low amplitude flutter that is seen on many 
subtraction images. 

An established technique to adaptively correct for both 
intra and inter-view motion is the navigator echo(^ 
acquisition. 29 In this method, an additional echo^NAVXls 
acquired per phase encoding echo. Since the TR is short 
(—10 ms) for the 3D SPGR sequence used for this study, 
acquisition of the additional navigator echos would substan¬ 
tially increase the scan time, thereby reducing the dynamic 
temporal sampling rate. Alternatively, a single navigator 
echo can be acquired each, along *y = 0 and kz=* 0, with 
kx—0 being acquired as part of the dynamic acquisition. 

This would add negligible time to the 3D acquisition, main- 
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gaining the desired high temporal sampling rate. In effect, 
this approach makes the assumption that a single set of dis¬ 
placements in x % y t and z can summarize the motion artifact 
drver the entire volume acquisition. This is similar to the 
basic assumptions described in this paper. Furthermore, since 
a displacement in space corresponds to a linear phase roll in 
&-space, the phase estimation technique described in this pa¬ 
per is considered appropriate to detect and correct the ex¬ 
pected motions. The averaging over multiple echoes in the 
dataset serves to provide some degree of noise insensitivity. 

In summary, the motion correction algorithm presented in 
this paper successfully reduces gross translation^ motion by 
estimating an averaged linear phase deviation per spatial fre¬ 
quency axis, between the reference and dynamic datasets. 
Where motion was negligible or could be attributed to other 
sources, implementation of the algorithm was not detrimen¬ 
tal to the original data. If the nature of motion was within the 
parameters of the model, fairly good correction in terms of 
improved visualization of structures of interest, reduction in 
subtraction edges, and suppression of background ghost arti¬ 
fact and blurring was observed. 
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